Exploring spin-orbital models with dipolar fermions in zig-zag optical lattices 
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Ultra-cold dipolar spinor fermions in zig-zag type optical lattices can mimic spin-orbital models 
relevant in solid-state systems, as transition-metal oxides with partially filled d-levels, with the 
interesting advantage of reviving the quantum nature of orbital fluctuations. We discuss two different 
physical systems in which these models may be simulated, showing that the interplay between lattice 
geometry and spin-orbital quantum dynamics produces a wealth of novel quantum phases. 



Orbital degrees of freedom of electrons play an im- 
portant role in the formation of various novel phases 
observed in transition-metal oxides with partially filled 
d- levels [3, 0]. In Mott insulators, they may enhance 
thermal as well as quantum fluctuations Q and lead 
to spin/orbital liquid states 043' an d to spontaneously 
dimerized states without any long-range magnetic or- 
der [1, Q. In many solid-state systems, the orbital 
dynamics is often quenched, due to the coupling of 
the orbitals to Jahn- Teller phonons, and the study of 
the quantum nature of orbitals demands systems with 
strong super-exchange coupling between spins and or- 
bitals. However, in real materials, the coupling strengths 
are fixed by nature, being very difficult to modify, thus 
limiting the experimental access to a potentially vast 
phase diagram. 

Ultra-cold spinor gases in optical lattices open new fas- 
cinating perspectives for the analysis of the quantum na- 
ture of orbitals T3-12|. The coupling constants in these 
systems can be easily controlled by modifying the lattice 
parameters and/or by means of Feshbach resonances [liij ]. 
Different lattice geometries, including frustrated lattices, 
like triangular [yj and Kagome lattices [l5|, may be 
created by combining different counter-propagating laser 
beams, and superlattice techniques. In addition, not 
only the physics in the lowest band but also that in 
higher bands may be controllably studied [l6j]. More- 
over, recent experiments on Chromium 17| and Dyspro- 
sium [lH atomic gases, and polar molecules [19 ]. are un- 
veiling the exciting physics of ultra-cold dipolar gases, for 
which the dipolc-dipole interactions may lead to exotic 
phases 

In this Letter, we show that dipolar spinor Fermi gases 
in appropriate zig-zag lattice geometries allow for the 
quantum simulation of spin-orbital models for a family 
of Mott insulating compounds, including systems with 
weak 



22|, |23[ as well as pronounced [24[ relativis- 
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tic spin-orbit couplings. Moreover, these models, which 
are relevant for real materials, as pyroxene titanium and 
layered vanadium oxides, may be explored with dipolar 
Fermi gases in parameter regimes which are hardly ac- 
cessible for solid-state compounds, allowing for the ob- 
servation of novel quantum phases. 

Physical realizations and effective Hamiltonian- We 




FIG. 1: (Color online) Two different systems discussed in this 
work, (a) Snake-like ladder with rungs (bonds) along y (x), 
and one fermion per rung. The inter (intra)-rung hopping is 
indicated by t (A); (b) Zig-zag lattice with one fermion per 
site occupying either p x or p y orbital states. Along a bond 
in x (y) direction only p x (p y ) orbitals are connected by a 
finite hopping amplitude t. A deformation of the lattice wells 
(shaded area) leads to a mixing between p x and p y orbitals 
with amplitude A. 



start by presenting two possible scenarios, A and B, in 
which dipolar Fermi gases may allow for the simulation 
of the above-mentioned spin-orbital models. 

A: First scenario is provided by a snake-like ladder 
as that shown in Fig. QJa), formed by bonds and rungs 
along x and y directions, respectively. This lattice geom- 
etry may be achieved by using a combined blue-detuned 
ladder-like lattice and a properly aligned red-detuned 
zig-zag lattice j25|. The overimposed red-detuned lat- 
tice allows for enhancing the potential barrier at alter- 
nated bonds in the upper and lower legs. A sufficiently 
strong red-detuned lattice leads hence to broken bonds 
as depicted in Fig. HJa). A similar technique has been 
recently employed for the realization of a Kagome lat- 
tice [15|. The hopping within the same rung is denoted 
by A, whereas the inter-rung hopping for unblocked bonds 
is denoted as t [Fig.QJa)]. The on-site repulsion U results 
from the combination of dipolar and contact interactions. 
The dipoles are oriented in the rrz-plane in such a way 
that fermions at the same rung experience maximal non- 
local repulsion V and dipoles on neighboring rungs inter- 
act identically, irrespective whether they occupy upper 
or lower sites, with the interaction strength much weaker 
than intra-rung repulsion. As a result, inter-rung dipolar 
interaction plays no role in the discussion below (inter- 
action between next-nearest rungs is considered negligi- 
bly small). We assume U, V » t, A, and consider the 
case of one fermion per rung. The system is then in the 
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Mott-insulator regime with one fcrmion localized on each 
rung. The sites belonging to the i-th rung of the ladder 
are distinguished by the pseudo-orbital quantum num- 
ber erf, with the convention that of = +(— )1 when the 
upper (lower) site on a given rung is occupied. Defin- 
ing a = U/V, and setting t 2 /2U as the energy unit we 
arrive to the effective spin-orbital Hamiltonian of Kugel- 
Khomskii type (26. 27|. 



H = ^[2S 2 S l+1 + a-i][l + (-l)Vf][l + (-l)Vr +1 ] 

i 

N N 

- A^2S i S i+1 [l- -?<r|+ 1 ]-A^ -f, (1) 



where Si is the spin-i operator (stemming from the 
spinor nature of the Fermi gas), and o~^ x are Pauli ma- 
trices describing the pseudo-orbital variables. We have 
added an additional term proportional to the coupling 
constant A. Although A = for realization A, it plays 
an important role in the alternative case B discussed be- 
low. Note that the ratio a may be modified basically 
at will, since U and V may be independently controlled 
using Feshbach resonances [13] , altering the lattice spac- 
ings, or modifying the transversal confinement (2p| . 

B: Second possible scenario is provided by a zig- 
zag lattice in the xy-plane, loaded with spinor dipolar 
fermions in p-bands [16( [Fig. [Tfb)], where the dipoles 
are oriented along z-axis. Assuming a strong confine- 
ment along z, we retain two degenerate orthogonal p x 
and p y orbitals per lattice site. In this realization, t de- 
notes the hopping between similar orbitals at neighboring 
wells [Fig. [Hb)]. An in-plane deformation of the lattice 
wells (e.g. by an additional weak tilted lattice) leads to 
a mixing of the p x and p y orbitals within the same well 
with an amplitude A. 

The interaction parameters for two fermions within the 
same well are on-site repulsions U (within the same or- 
bitals) and V (among different orbitals), and Hund's ex- 
change Jh [27| . Two fermions occupying the same or- 
bital may form a symmetric or an antisymmetric state 
with respect to the orbital index with corresponding en- 
ergies U+Jr and U — Jh, which are split by the so-called 
pair- hopping term with an amplitude Jh |27j ■ The over- 
all energy scale will now be modified to t 2 /2U, where 
U = (U 2 — J H )/U. When two fermions occupy dif- 
ferent orbitals, they may form a spin-singlet or a spin- 
triplet state with corresponding energies V + Jh and 
V — Jh, which are split by Hund's exchange. In the 
strong coupling limit U ± Jh, V ± Jh ^ t, A and with 
one particle per lattice well the system is in the Mott- 
insulator regime, and we arrive at Hamiltonian ([T]) with 
A = J H U/(V 2 - Jh) [27j]. Now, a also gets modified 
accordingly, a = U(V + J H /2)/(V 2 - J 2 H ). 

Note that for a purely contact interaction V — Jh and 
in the triplet channel two fermions do not experience any 



repulsion [27j. Thus, without dipolar couplings, the Mott 
phase of one fermion per well would not be stable. In the 
Mott-insulator regime, one can still vary A in a wide 
range by changing the relative ratio of the strengths of 
the contact and dipolar interactions. 

The model (JTJ, at a ~ 1, describes the spin-orbital in- 
terplay in Mott insulating transition-metal compounds. 
At A = 0, it describes a zig-zag chain of spin one-half 



Ti 3H 



ions, with active d xy and d yz orbitals, in pyroxene 



titanium oxides ^TiSi 2 6 (A = Na,Li) |22j, |23|. On the 
other hand, for A = the Hamiltonian ([1]) represents 
the ID counterpart of the 2D model for Sr2V04 24 [. 
In the latter, the role of spins in Eq. (QJ is played 
by an isospin variable discerning the Kramers partners, 
while the pseudo-orbital variables distinguishes two low- 
est Kramers doublets of V 4+ ion, and A term in Eq. ([T]) 
represents relativistic spin-orbit coupling. 

In the following we study the ground-state properties 
of the model ([1]) in the wide parameter range by analyt- 
ical and complementary numerical approaches. 

Ground-state phase diagram for A = 0.- We first dis- 
cuss the case A = relevant for realization A. At A = 
the orbitals become classical and the ground-state phase 
diagram is easily mapped. For a > 2, the ground- 
state is 2 x 2 N degenerate: there is anti-ferro (AF) or- 
der in orbitals (erf) = ±(—1)' (± refers to two degen- 
erate AF states), whereas the spin part is completely 
degenerate (bold line on Fig. [UJa)). For a > 2 and 
A — >• 0, the system is described by an effective spin model, 
Hs ~ A z2i SiSj+i, an isotropic Heisenberg antifcrro- 
magnet (iH). Thus, quantum fluctuations of orbitals, in- 
duced by A-term, lifts immediately the infinite degen- 
eracy of the ground state, resembling order from dis- 
order. We denote this phase (iH,AF), where the first 
term denotes iH spin phase, and the second term AF 
orbital phase. We employ a similar notation from now 
on. Higher order terms in A/a cannot break the SU(2) 
spin and translational symmetries, and thus for A — > 
the iH phase in spin degrees is stable. As shown be- 
low, the iH is recovered for strong A independently of the 
value of a. Thus, we can expect that for a > 2 there 
is an unique iH phase in spin degrees of freedom for any 
A 7^ 0. On the contrary, with increasing A the orbital de- 
grees of freedom experience an Ising transition from AF 
to paramagnetic (P) phase with (erf ) = 0. 

For A = and < a < 2, the exact ground state 
is two-fold degenerate and represents a direct product 
of ferro (F) orbital order, (erf) = +1 (—1), and spon- 
taneously dimerized Majumdar- Ghosh (MG) state [111 
in spins, with spin-singlets located on odd (even) bonds. 
We call this phase a dimer- ferro (D,F). An infinitesimal A 
generates an exchange between the disconnected nearest- 
neighbour dimers, ~ A 4 Y^, ^2i+i^2i+2- With increasing 
A, the dimerization order in spins, D = J^i KS<S»+i — 
SjSj_i)|, disappears together with the orbital ferro or- 
der at a Kosterlitz-Thouless (KT) phase transition (29| . 
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FIG. 2: Ground state phase diagram of the model JTJ) for (a) 
A = and (b) a — 1. First and second phases in parentheses 
refer to spin and orbital sectors, respectively. See text, for a 
description of phases and phase transitions. 



A numerical ground-state phase diagram of the model ([T]) 
for A = is presented in Fig. (JJa) (details of the numer- 
ical simulations are discussed later). 

Ground-state phase diagram for A ^ 0.- We now turn 
to the case of finite A. We focus on the regime a = 1 rele- 
vant to realistic condensed- matter systems [1, 22 -24 1. For 
A = 0, a simple calculation shows that the ground state 
is (F,AF) for A > 2(2 - a) and (D,F) for A < 2(2 - a). 
Note, that once the spin dimerization pattern is sponta- 
neously chosen, the direction of orbitals becomes unam- 
biguously selected, thus Ising Z 2 orbital order is 'slaved' 
by translation symmetry breaking. Whereas for A = an 



infinitesimal A induces AF spin-exchange oc A on inter- 
dimer (weaker) bonds in the (D,F) phase, for finite A 
the leading spin-exchange along the weaker bonds is in- 
stead ferromagnetic ~ — X 2 AJ2 82^+182^+2- Hence there 
is a competition between the A and A terms promot- 
ing, respectively, a F and AF character of the weaker 
bonds. The character of spin correlations on weak bonds 
(821+182^+2) changes from AF (small A region) to F 
(larger A region) across Amg = A 2 /4 + 0(A 4 ) line, re- 
sembling the behaviour across the MG point in the spin-i 
h ~h model 0. 

For A —> 00 the orbital degrees of freedom are 
quenched. The system becomes equivalent to a SU(2)- 
symmetric spin-i chain, which up to 0(\~ 2 ) is described 
by the effective Hamiltonian, 

h s = £ [inSA+n + r!(s l _ 1 s i )(s 4+1 s J+2 )] (2) 



where ji = 2 - 2A + 0(A _1 ), j 2 = A -1 , j 3 , n are both of 
order A -2 , and the longer range exchanges are suppressed 
as j n ~ 0(A 1_n ) [27]]. There are two clear phases in this 
regime, (iH,P) for 1 - A > A" 1 and (F,P) for I — A <C 
-A" 1 . 

One may suspect for large A a direct (iH,P) to (F,P) 
transition with growing A. This, however, is not the 
case, as can be shown by a bosonization analysis of the 
effective spin model ©. Starting from (iH,P) state, with 



increasing A the system necessarily enters first into a 
dimerized state via a KT phase transition at Akt — 
1 — 1/2A. At larger A, a MG state will be approximated 
at Ajv/g — 1 + 1/2A, where dimerization will reach the 
value D ~ 3/4. Increasing further A, at A' ~ 1 + 3/2A, 
nearest-neighbor coupling in Eq. $2$ vanishes, ji(A') = 0, 
and the 1/A 2 terms, 7 3 an d f2, are the leading ones that 
couple two sub-chains |27l| . Bosonization shows [27| that, 
inspite of these terms, in the large A limit, the system be- 
haves at low energies as two decoupled spin-^ chains, al- 
beit at A* = A'+(9(A" 2 ). Hence around the A = A* line 
the effective spin model in strong coupling is described 
by a ji — J2 model, where j\ changes sign from antifcrro 
(for A < A*) to ferro (for A > A*), whereas j'2 stays 
positive. Bosonization, supported by recent numerical 
studies, predicts that the ground state of two weakly cou- 
pled spin-i chains is dimerized irrespective of the sign of 
ji coupling 3lM33| . Hence, in our case there is a special 
fine-tuning line bisecting the dimerized phase, A = A* 



line, described by double KT phase transition [34j where 
spin dimerization and ferro orbital order both vanish. Fi- 
nally, there is a first order phase transition line separating 
(D,F) and (F,P) states at A F ~ 1 + 7/2A. 

This sequence of phases and phase transition curves 
(Akt, A* , Af), which has been established analytically 
for large A with the help of effective spin model ©, has 
been confirmed for A > 4 by numerical simulations (dis- 
cussed below) of the original model ([1} . 

Numerical procedures - The phase diagrams depicted 
in Figs. [5Ja,b) were obtained by means of a combination 
of Lanczos exact diagonalization, density matrix renor- 
malization group simulations based on matrix product 
states (MPS) [35J and the infinite time-evolving block 
decimation (iTEBD) algorithm j3f|, confirming the ana- 
lytical predictions for A -+ and A -+ 00. 

The KT transition between (D,F) and (iH,P) was ex- 
tracted by Lanczos method from the extrapolation of the 
level crossing [37[ between the first excited singlet of the 
(D,F) phase and the first excited triplet of the (iH,P) 
phase for systems of up to N = 12 rungs/ wells. Ising 
transition lines were obtained, by MPS simulations, from 
the peak in the fidelity susceptibility [38|, and that on 
Fig. [2Jb) between (F,AF) and (F,P) accurately follows 
the analytical line A = A/2 + a. The first order phase 
transition into the (F,AF) or (F,P) states [bold line in 
Fig. [2](b)] was obtained by MPS method, from the jump 
of the ground-state total spin from (singlet state) to a 
fully polarized N/2 (ferro state) [3^. In MPS simulations 
we have used periodic boundary conditions and system 
sizes of up to N = 32 rungs/wells were considered. For 
A > 4, the double KT phase transition line bisecting the 
(D,F) phase [see Fig. [2][b)] was determined from vanish- 
ing dimer order using iTEBD method. For smaller values 
of A (1 < A < 4), around the dashed line of Fig. [1(b), 
our simulations indicate an intermediate quadrumerized 
phase 0. 
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Final remarks - We have assumed above a unit occu- 
pation per rung/well, for which dipolar interactions were 
necessary for realizing Mott insulating state. The case of 
two fermions per rung/well does not require dipolar inter- 
actions. In the case of the snake-like lattice [see Fig.QJa)] 
with two fermions per rung, a dimerized state along the 
rungs (for A > t) will be separated from a dimerized 
state along the bonds (A < t) by a KT phase transition 
at A = t. If in the zig-zag geometry of Fig. [ljb) we place 
two fermions per well, both mp orbitals, then for A <C Jh 
due to the Hund's coupling a total S = 1 state will be 
formed in each well, and the Haldane phase of a spin-1 
chain will be realized. 

In conclusion, dipolar fermions on zig-zag lattices 
can capture relevant spin-orbital models of realistic d- 
electron systems and allow to explore parameter regimes 
which are hardly accessible for solid-state compounds. 
Moreover, the quantum nature of orbital fluctuations can 
be revived, which combined with geometric frustration 
and spin dynamics produces an intriguing rich ground 
state phase diagram. 
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SUPPLEMENTARY MATERIAL TO 
"EXPLORING SPIN-ORBITAL MODELS WITH 
DIPOLAR FERMIONS IN ZIG-ZAG OPTICAL 
LATTICES" 

In this supplementary material, we provide addi- 
tional details concerning the derivation of the spin-orbital 
Hamiltonian and the bosonization procedure. 



THE EFFECTIVE SPIN-ORBITAL 
HAMILTONIANS 

In this section we outline the essential steps of the 
derivation of the effective spin-orbital Hamiltonian for 
the two systems, A and B, discussed in the Letter. 



A: Snake-like lattice 

We introduce the fermionic annihilation operator c a i jS , 
where a = 1 (2) indicates the up (down) sites on i-th 
rung, and s =t or I refer to the spin. Using the notation 
of the Letter, system A is described by the Hubbard like 
Hamiltonian H = Hk in + H int , where 



(3) 



i.a.s 

accounts for single-particle processes, and 

H in t = UY^cii,^ai,^ aU c aiA (4) 

2, a 

+ V C li,s C U <S ct itS ,C2i, s >. 

corresponds to two-particle interactions. In the limit 
[/, V 3> t, and retaining one particle per ladder rung 
(quarter filling), the effective Hamiltonian up to the sec- 
ond order in t acquires the form of Eq. (1) of the Let- 
ter (at A = 0), where 

a a 

Sf = ^ = ^ c ai^ c ai,\. i 
a a 

are the spin-i operators on the i-th rung, and 

S 

s 

are the Pauli matrices corresponding to the pseudo- 
orbital degrees of freedom. 



B: p-band zig-zag lattice 

In the system B, the pseudo orbital index a = 1 (2) 
refers to p x (p y ) orbital in a given well. The kinetic part 
of the Hamiltonian remains the same as in Eq. ([3]), while 
the interaction part, in addition to the terms shown in 
Eq. ((4]), acquires the following supplementary terms: 



n int 



(5) 



J H Yl I c i «T c i 4 C2 *t + 4 *l4 4 Cl U Cl <t] . 



where Si( 2 )i and rii(2), are, respectively, spin and den- 
sity operators for a fermion in p x (j> y ) orbital state. The 
first line stands for Hund's exchange and the second one 
describes the so-called pair- hopping 

The coupling constants U, V, and Jh, entering in 
Eqs. (|4I5[) . can be expressed in terms of a two-body po- 
tential and single-particle wave functions One finds: 



U = 
V 



J dr 1 dr 2 p2 (y) (r 1 )U(r 1 - r 2 )P^ y) (r 2 ), (6) 
J dnd^P^n)!/^ - r 2 )P y 2 (r 2 ), (7) 



and 

J H = /dndraP^rOP^n^n - r 2 )P,(r 2 )P y (r 2 ). (8) 

Above P x and P y are orbital wavefunctions of the same 
well, and V{y\ — r 2 ) is the total intcrparticle potential 
(including both contact as well as dipolar interactions). 

The energy of the state with two fermions occupy- 
ing the same orbital is U . However, the pair-hopping 
term moves two particles from the doubly occupied or- 
bital to an empty orbital with the amplitude Jh, and 
true eigenstates become orbital-symmetric and orbital- 
antisymmetric states with corresponding energies U + Jh 
and U — Jh- 

Two fermions occupying different orbitals may form 
a spin-singlet or a spin-triplet state with corresponding 
energies V + Jh and V — Jh 0] ■ 

In the limit U ± Jh, V ± Jh 3> t, and up to second 
order in t, the effective Hamiltonian acquires the form 
H = J^i Hi : i+\ : where P^i+i is the two-site Hamiltonian: 



Hii+l 



t 



^P M+1 (S T =0)(f + (-l)V*)(l + (-I)V +1 ) 
P M+1 (S T =0)(l-<a,f +1 ) 



2(V+ J H ) 
t 2 



2(v _ jH f^(s T =m-^ +1 ) 



(9) 
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FIG. 3: Two sub-chains with intrachain coupling j% (bold 
lines) coupled by zig-zag interchain coupling j\ (dashed line). 

where we have introduced the operators Pi i i+i(S T =0) = 
-SiSi+i + 1/4, and P %l+1 (S T =l) = S^+i +3/4, which 
project onto two-fermion states on sites i and i + 1 with, 
respectively, total spin S T = and S T = 1, and U = 
(U 2 - J%)/U. 

The first line in Eq. ([9]) accounts for the situation 
in which both fermions occupy same orbitals on sites 
i and i + 1. In that case, when a fcrmion hops to the 
neighboring site Q, the orbital-symmetric or orbital- 
antisymmetric states of the spin-singlet pair is reached 
(note that the formation of the spin-triplet pair is for- 
bidden by Pauli exclusion). 

The second and third lines of Eq. (0) account for the 
configuration in which two fermions occupy orthogonal 
orbitals on neighboring sites. When a fcrmion hops to 
a (singly) occupied neighboring site (remembering that 
inter-site hopping t does not change orbital quantum 
number) , depending on the total spin of the two fermions 
the interaction may take place either in the spin-singlet 
channel, with energy cost V + Jh, or in the spin-triplet 
channel, with energy cost V — Jh- A proper re-ordering 
of the terms in Eq. © leads to the Hamiltonian of Eq. (1) 
of the Letter. 



BOSONIZATION ANALYSIS 

In the following we provide some details on our 
bosonization analysis of the large A regime (A — > 00), and 
in particular on how this approach provides an under- 
standing of the special fine-tuning line A = A* cutting 
the (D,F) phase at which spin dimerization and orbital 
ferro orders vanish. 

Integrating out the orbitals in large A limit, to the 
second order in 1/A we obtain the effective spin-i model 
given in Eq. (2) of the Letter with, 

J1 = 2 (l-A) + 4 + A ^ + A ) + ( A^ 

h = \ + 0(X-% , 3 = i±^ + 0(A- 3 ), 

Q = ^ + 0(\- 3 ) , (10) 

where we use the same notations as in the Letter. A — > 00 
is a convenient limit for bosonization, in particular for 
ji *C ,]2- In that case, one can consider two weakly cou- 
pled chains (see Fig. [3]), which we denote as Sj = S21 and 



Sf = S2i+i. Using bosonization identification, the spins 
along each chain are represented by smooth and stag- 
gered parts Sf -> 3% + 3% + (-iyn a Q. The coupling 
of the two chains induces marginal operators. Apart 
from the scalar (dimerization) operators, there are spin- 
full marginal terms, so-called twist operators @. Up to 
irrelevant and marginal couplings that induce only ve- 
locity renormalization, the effective bosonized Hamilto- 
nian density describing the spin degrees of freedom in the 
large- A limit is of the form 



n = — J|J| + A(J[J 2 fl + jyi) + D 2 3 a L 3 a R (11) 

+T 1 e ab (e a d x e b + n a d x n b ) + T 2 £ ab (3e a d x e b - n a d x n b ). 



where v ~ j'2, a summation convention is assumed for 
repeated indices a,b — (1,2), and £ = (L, R), e a b is the 
antisymmetric symbol, and e a stand for the dimerization 
operators, such that S"S° +1 ~ ( — l) 1 e a (x)+ a less rele- 
vant smooth part. 

Bare values of the twist couplings 7i 2, and the inter- 
chain dimerization D\ depend on and O, in linear 
order (with non- universal proportionality coefficients), 
whereas the intra-chain dimerization amplitude D2 de- 
pends on J2 and f2. 

One-loop RG equations corresponding to the effective 
model pT|) are identical to those obtained in Q for the 
SU (2) symmetric case. RG flow is dominated by inter- 
chain dimerization D± , which in the infrared limit always 
scales to strong coupling, except for the initial condition 
T\ = T2, corresponding to the fine tuning A = A*, which 
we interpret as an infrared decoupling point of the two 
chains. In contrast D2 scales to zero at low energies for 
A — > 00. Thus for A — > 00 there is a special line in 
(D,F) phase where spin dimerization (and consequently 
the ferro orbital order) vanishes. 

In fact, bosonization is only needed to capture the in- 
fluence of j3 and f2 terms close to j\ = 0. Away of 
that region, for A —5- 00, one just needs to retain ji and 
j'2 terms, and borrow known results from the frustrated 
ji — ]2 spin-i chain. This gives us an estimate of other 
phase transition lines (we use the notations of the Let- 
ter) Akt and A p which follow respectively from ji ~ ij 2 
and ji = — 4j2 conditions. Whereas Amg (which does 
not represent a phase transition line) follows from the 
condition j\ = 2j2- 



* Also at Andronikashvili Institute of Physics, 0177 Tbilisi, 
Georgia. 
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